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Abstract 


An analysis of the non-modal growth of perturbations in a boundary layer in the 
presence of a streamwise pressure gradient is presented. The analysis is based on PSE 
equations for an incompressible fluid. Examples with Falkner-Skan profiles indicate that 
a favorable pressure gradient decreases the non-modal growth while an unfavorable 
pressure gradient leads to an increase of the amplification. It is suggested that the 
transient growth mechanism be utilized to choose optimal parameters of tripping 
elements on a low-pressure turbine (LPT) airfoil. As an example, a boundary-layer flow 
with a streamwise pressure gradient corresponding to the pressure distribution over a LPT 
airfoil is considered. It is shown that there is an optimal spacing of the tripping elements 
and that the transient growth effect depends on the starting point. At very low Reynolds 
numbers, there is a possibility to enhance the transient energy growth by means of wall 
cooling. 


Nomenclature 


E 

G 

H 

H l 

H u 

H re f 

L 

Lref 

L s 

Rq 

Re/, 

T w / T ad 

U 


= energy norm 
= E out IE in energy ratio 
= shape factor 
= JvL/U eL 

~ / G ex it 

= scale in v- and z- directions 
= length along the streamwise direction 
= reference length 
= surface length 
= Reynolds number 
= U eL L / v 

= temperature factor (ratio of the wall temperature to the temperature of the 
adiabatic wall) 

= streamwise velocity component of the mean-flow velocity 
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U eL 

Uref 
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= free-stream velocity at x = L 
= reference velocity 

= nomial-to-the wall velocity component of the mean-flow velocity 
= streamwise velocity disturbance 
= nonnal velocity disturbance 
= spanwise velocity disturbance 
= streamwise coordinate 
= coordinate nonnal to the wall 
= spanwise coordinate 
= spanwise wave number 
= Hartree parameter 

00 ^ ^ ref L ref 

= kinematic viscosity 
= momentum thickness 
= pressure disturbance 


Superscripts 
T = transposed 


Subscripts 

e = free stream 

exit = exit conditions 

in = starting point 

out = ending point 


1. Introduction 

Laminar-turbulent transition in shear flows is still an enigma in the area of fluid 
mechanics. The conventional explanation of the phenomenon is based on the instability 
of the shear flow with respect to infinitesimal disturbances. The conventional 
hydrodynamic stability theory deals with the analysis of normal modes that might be 
unstable. The latter circumstance is accompanied by an exponential growth of the 
disturbances that might lead to laminar-turbulent transition. Nevertheless, in many cases, 
the transition scenario bypasses the stage of the exponential growth associated with the 
nonnal modes. This type of transition is called as bypass transition. Observations of 
laminar-turbulent transition in plane Couette flow, in circular pipe flow, in boundary 
layers at relatively low Reynolds numbers, etc., can serve as examples of bypass 
transition. An understanding of the phenomenon has eluded us to this day. One 
possibility is that bypass transition is associated with so-called algebraic growth of 
disturbances in shear flows . 1 ” 3 

The phenomenon of algebraic growth of disturbances in shear flows has been of great 
interest during the last two decades as it may be associated with the bypass transition 
mechanism. Ellingsen and Palm 1 considered, in the inviscid case, an initial disturbance 
independent of the streamwise coordinate and found that the streamwise disturbance 
amplitude may grow in time, even though the basic flow does not posses an inflection 
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point. Landahl" showed that all parallel inviscid shear flows are unstable to a wide class 
of three-dimensional disturbances. The result is independent of whether or not the shear 
flow is unstable to an exponential growth of wavelike disturbances. This type of 
instability that is not related to exponential growth is also referred to as “ non-modal 
growth. ” Mathematically, the effect of non-modal growth is associated with non- 
normality of the linearized Navier-Stokes operator and non-orthogonality of the 
eigenfunctions. In simple words, the essence of non-modal growth is the possibility of 
combining the exponentially decaying modes in such way that their sum will possess 
transient growth. One can find a vast bibliography on the topic in the Otto Laporte Award 
Lecture by Reshotko 3 and in a monograph by Schmid and Henningson. 4 

Numerical analysis of spatial non-modal growth within the scope of the linearized 
boundary-layer equations for an incompressible flow over a flat plate was carried out in 
Refs. 5 and 6. Spatial analysis within the scope of the linearized Navier-Stokes equations 
(quasi-parallel approximation of compressible and incompressible flows) was presented 
in Refs. 7 to 9. The main results of these theoretical models are as follows: 

• A system of counter-rotating streamwise vortices, which are periodic in the 
spanwise direction, provides the strongest growth of the disturbance. 

• There is an optimal spacing of the streamwise vortices, leading to the strongest 
effect. 

The effect of pressure gradients on the transient growth mechanism was considered 
within the scope of temporal theory by Corbett and Bottaro 10 and within the scope of 
spatial theory by Tumin and Reshotko. 9 These results were based on the quasi-parallel 
flow assumption. Turnin’ 1 analyzed the pressure gradient effect for the Falkner-Skan 
profile within the scope of an analytical model when the spanwise wave number is very 
small. The pressure-gradient effect within the scope of spatial theory with nonparallel 
base flow and finite spanwise wave numbers has not been considered, yet. 

Another motivation for the present work stems from separation flow control on low- 
pressure turbines (LPT). The perfonnance of LPTs is strongly affected by the flow 
separation. There is a possibility of delaying the boundary-layer separation by tripping 
the boundary layer with the help of roughness elements or other devices. Usually, trial- 
and-error method is used to determine an appropriate placement of the control elements. 
This approach is time consuming and expensive. A recent investigation by Reshotko and 
Turnin’ 2 demonstrated that roughness-induced transition might be related to the transient 
growth mechanism. 

Periodically spaced in the spanwise direction, roughness elements generate a system 
of counter-rotating streamwise vortices. Due to a secondary instability mechanism, the 
streamwise vortices can lead to earlier transition to turbulence. They also provide a 
mixing enhancement due to redistribution of the streamwise momentum. Consequently, 
optimization of the streamwise vortices for maximum energy growth leads to maximizing 
of the flow control effectiveness. In the present work, analysis of the optimal 
disturbances/streamwise vortices associated with the transient growth mechanism will be 
perfonned for boundary layers in the presence of a streamwise pressure gradient. The 
theory will provide the optimal spacing of the control elements in the spanwise direction 
and their placement in the streamwise direction. 
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2. Governing Equations 


In the present work, we utilize the method developed in Ref. 19. Because the flows of 
interest have relatively low Mach numbers, we consider steady three-dimensional 
disturbances in an incompressible two-dimensional boundary layer. We choose the 
streamwise coordinate x along the surface. The coordinate y will measure distance from 

the wall. We define a small parameter s = v/U f L 7, where v , U f , and L f are 

viscosity, reference velocity, and reference length, respectively. The streamwise 
coordinate is scaled with L ^ while the vertical coordinate y and spanwise coordinate 

z are scaled with ^j v ^ re f / U re j ■ The following scaling is assumed for the velocity 
disturbances u, v, and w , and the pressure n : 


u ~U 


ref ’ 


v ~ sU 


ref’ 


W~ sU 


ref ' 


P ~ £~pU t 


ref 


( 1 ) 


This scaling of the linearized Navier-Stokes equations and neglecting of the curvature 
effects lead to the governing equations for Gortler instability with the Gortler number 
equal to zero. We look for a periodic solution in the spanwise direction with the 
corresponding wave number J3. The governing equations for the amplitude functions can 
be written in dimensionless form as follows: 5 ' 6 


dudv 

1 1 - Bw= 0 

dx dv 


d ( s du dU d z u 2 

—(Uu) + V — + v— = — y~P U 
dx dy dy dy 


^(uV + vU) + ^-{2Vv) + flVw+^ = ^f-/} 2 v 
dx dy dy Qy 2 


2 

d t TT \ d / T r \ n d w n 2 
— ( Uw) + — (Vw)-/3p = — --p w 

dx v dy dy 2 


( 2 ) 

( 3 ) 

( 4 ) 

( 5 ) 


where U (x, y) and V (x,v) are the streamwise and normal velocity components of the 

base flow, respectively (in addition, the latter is divided by s). The following boundary 
conditions are applied to the solutions: 

v = 0: u = v = w = 0 (6a) 


y— >oo; u,w,?r 0 (6b) 

The equations (2) to (5) can be solved subject to boundary conditions (6a) and (6b) with 
prescribed initial velocity perturbations at x = xq. 
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The governing equations can be recast as follows: 

(Af)=B 0 f + B,f,.+B,f„ (7) 

where A, B Q , B| , and B 0 are 4x4 matrices (one can find them in the Appendix; see 

also Ref. 5) and f = (u,v,w,p) T . The superscript “T” stands for “transposed,” and the 
subscripts “x” and “v” denote differentiation with respect to x and y, respectively . 

3. Optimization of Energy Growth 

The authors of Refs. 5 and 6 employed an iterative procedure to find the optimal 
disturbances in tenns of the maximum of the energy growth ratio G = E out / E in , where 
Ej n and E out stand for the input and output energy nonns. Andersson et al. 5 used the 
same definitions of E in and E out as for the disturbance energy 



0 


whereas Luchini 6 employed the knowledge that the optimal disturbances are represented 
by streamwise vortices with corresponding output as the streamwise velocity streaks, 


p. 
n in 



ymax 

Eout = J u2 dy 
0 


G = 




y max 

| u 2 dy 
0 


| {v 2 +w 2 ^dy 
0 


(9a) 


(9b) 


(9c) 


As was shown in Ref. 5, the two definitions of the optimal disturbances lead to the 
same results at Reynolds numbers of 10 4 and higher. Because the iteration procedure 
based on the optimization of ratio (9c) provides significant simplification, we adopt it for 

9 

the following analysis. Because Eqs. (2) to (5) are independent of s , the value of s G is 
invariant with respect to the Reynolds number. 
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4. Numerical Results 


4.1 Falkner-Skan Base Flow 

We consider a Falkner-Skan family of boundary-layer profiles with free-stream velocity 
distribution U e = Cx m and corresponding Hartree parameter = 2m / ( m + 1) . For the 

purpose of convenience, we have used the velocity scale U re j = U eL = CL m and the 
length scale L re f = LI ( in + 1 ) . The latter allowed the use of the conventional scaling of 
boundary-layer solutions with H re j- = ^JvL / (in + \ )U eL = JvL re j- / U re j- . 

Figure 1 shows the scaled the energy ratio versus spanwise wave number (5 for three 
Hartree parameters, (3 H = -0.1, 0.0, and 0.1. The starting and the ending points, x in / L 
and x out / L , are equal to 0.2 and 1.0, respectively. The Reynolds number Re/ in Fig. 1 
and what follows is defined as U eL L / v . One can see that an unfavorable pressure 
gradient < 0) leads to an increase in the energy growth while a favorable pressure 

gradient (ftn > 0) leads to suppression of the transient growth mechanism. The latter is 
consistent with results obtained within the scope of parallel flow approximation. 3 



Fig. 1. Effects of the spanwise wave number [3 and the Hartree parameter on 
transient growth (starting point x in / L = 0.2 ). 
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Figures 2 and 3 show similar results, but the starting points are x in / L =0.4 and 0.6, 
respectively, respectively. A comparison of Figs. 1 to 3 indicates that there is a spanwise 
wave number, (3H re f = 0.5 -0.6, and a starting point, x in / L , that maximize the energy 

growth. These parameters correspond to optimal spanwise spacing and streamwise 
placing of perturbators for maximum flow control effectiveness. Figure 4 demonstrates 
the energy ratio versus the downstream coordinate xlL at fill re j =0.5 and x in /L = 

0 . 2 . 
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Fig. 4. Optimal energy ratio versus the downstream coordinate at three Hartree 
parameters ( x in I L = 0.2 , (3H L = 0.5 ). 
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4.2 Example No. 1 of LPT Conditions 

1 3 

Volino simulated low-pressure turbine (LPT) airfoil conditions in a low-speed wind 
tunnel. The test section was designed as a passage between two airfoils. The local free- 
stream velocity at a favorable pressure-gradient region was closely approximated by the 
following equation: 


Ue 

^ exit 


1.48 


f \0.214 
X 


vAs J 


( 10 ) 


where L s is the suction surface length and U exjt is the nominal exit free-stream velocity 
based on the inviscid solution. The distribution (10) corresponds to the Falkner-Skan 
flow with the Hartree parameter = 0.353 . 



Fig. 5. Effects of the spanwise wave number /? and the starting point x in I L on 
transient growth at conditions of the experiment in Ref. 13 (/?// = 0.353 , 

x out / L S 0*444). 

Figure 5 demonstrates the energy ratio scaled with the Reynolds number R e exit = 
U exit L s / v versus the spanwise wave number scaled with H Ls = JvL s / U exjt . The 
ending point was prescribed at x out / L s = 0.444 while the starting points varied from 
0.111 to 0.289. The streamwise velocity perturbation at x out / L s = 0.444 , 
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x in /L s =0.111, and J3H Ls = 0.925 is shown in Fig. 6, and the corresponding optimal 

profiles of v and w are presented in Fig. 7. 

The results indicate that we are dealing with a very strong favorable pressure gradient 
that suppresses the transient growth mechanism. If the Reynolds number is about 50,000, 
the transient growth will provide an energy amplification of less that 50. If we take into 
account that the real perturber will not produce the optimal inflow field, there is little 
likelihood that the bypass transition via a secondary instability of the velocity streaks will 
prevent the flow separation. 



Fig. 6. The streamwise velocity perturbation at the ending point x out / L s = 0.444. 
The parameters correspond to the experimental conditions in Ref. 13 ( x jn / L s = 

0.111, /3H Ls = 0.925). 
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Fig. 7. The optimal velocity perturbations at x in / L s =0.111 and j3H Ls = 0.925 
corresponding to the streamwise velocity perturbation at x out / L s = 0.444 shown in 
Fig. 6. 


There is a possibility of enhancing the transient growth mechanism by means of wall 
cooling. The effect of wall cooling was investigated by Tumin and Reshotko 9 within the 
scope of a parallel flow approximation. In order to estimate possible increases of the 
energy ratio on a cold wall at a high favorable pressure gradient, we utilize the method of 
Ref. 9 for a compressible flow with local Mach number 0.5 and Hartree parameter 0.353. 
The results are shown in Fig. 8. One can see that cooling of the wall might provide a 
tenfold increase in the energy ratio. 
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Fig. 8. Effect of the temperature factor on energy growth at the experimental 

conditions of Ref. 13. 


4.3 Example No. 2 of LPT Conditions 

Experiments on a flat plate with a pressure gradient corresponding to a modern low- 
pressure airfoil were carried out by Volino and Hultgren 14 at the NASA Glenn Research 
Center. Figure 9 shows a comparison of the theoretical and experimental free-stream 
velocity distributions along the flat plate (one can find more details in Ref. 14). Although 
there is a discrepancy between the theoretical and experimental data for the region of the 
favorable pressure gradient, we utilize a theoretical distribution in the following 
evaluations. The numerical methods for solving the boundary layer equations are 
described elsewhere. 15,16 Figures 10 and 11 show comparisons of the theoretical and 
experimental shape factors and the momentum thickness, respectively. One can see that 
the theoretical distribution of the free-stream velocity provides reasonable agreement 
between the theoretical and experimental data for the boundary-layer parameters. 
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Fig. 9. Comparison of theoretical and experimental 14 free-stream velocity 

distribution along the plate. 



Fig. 10. Comparison of theoretical shape factors at two Reynolds numbers with the 

experimental data in Ref. 14. 
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Fig. 11. Comparison of the theoretical momentum thickness with the experimental 

data in Ref. 14. 



Fig. 12. Comparison of the theoretical and experimental 14 velocity profiles 

at x / L re j = 0.266. 


Comparisons of the theoretical and experimental velocity profiles are shown in Figs. 12 
and 13. Because there are some discrepancies between the theoretical and experimental 
data, one should consider the following transient growth analysis as an estimation under 
the experimental conditions. 
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Fig. 13. Comparison of the theoretical and experimental 14 velocity profiles 

at x/ L re j = 0.51. 



Fig. 14. Effects of the spanwise wave number, /?, and the starting point, 
x in / L re f , on the energy growth at the experimental conditions of Ref. 14. 


Figure 14 demonstrates the energy ratio versus the spanwise wave number at 
x out / As =0.5 and three starting points, x in / L s = 0.125, 0.250 and 0.375. One can see 
that the results are similar to the previous case with a Hartree parameter of 0.353. 
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5. Summary 


The results for the transient growth phenomenon within the scope of the linearized 
boundary-layer equations in the presence of a streamwise pressure gradient are consistent 
with previous results obtained within the scope of the parallel flow approximation and 
linearized Navier-Stokes equations. 9 A favorable pressure gradient decreases the non- 
modal growth while an unfavorable pressure gradient leads to an increase of the 
amplification. 

The example of a Falkner-Skan flow with a Hartree parameter /3 H = 0.353 
corresponds to the experimental data and simulates the flow over a low-pressure turbine 
airfoil upstream of the separation point. At this pressure gradient, the transient growth 
mechanism appears to be weak, and the latter, together with the low Reynolds number 
allows us to speculate that tripping of the boundary layer via a secondary instability of 
streaky structures is likely to fail. This conclusion should be verified by further work. 
Estimation of the transient growth at experimental conditions on a flat plate with a 
streamwise pressure gradient 14 shows that the results are similar to the example of 
Falkner-Skan flow at /3 H = 0.353 . The theory of the transient growth mechanism 

predicts the possibility of enhancing the energy growth by means of wall cooling. The 
example within the scope of the parallel flow theory 9 demonstrates that cooling of the 
wall might provide a tenfold increase in the energy ratio. Future experiments on boundary 
layer tripping accompanied by wall cooling will contribute to our understanding of the 
bypass transition mechanism. 

Theoretical analysis of the optimal disturbances predicts that there is an optimal 
spacing between perturbers and their optimal location from the leading edge. The latter 
results can be utilized in future experiments with tripping of the boundary layer over the 
LPT airfoil. This type of experiment can be carried out at the NASA Glenn Research 
Center, 14 at the US Naval Academy, 13 and at the University of Notre Dame. 17 

Consideration of the optimal velocity perturbations in Fig. 7 indicates that they are 
spreading across the boundary layer. This means that an array of generators localized on 
the wall will not provide excitation of the optimal disturbances. Therefore, the question 
of realizability of the optimal disturbances arises. For example, one can solve the 
receptivity problem for an array of generators on the wall and evaluate their shapes (or 
other parameters), providing velocity disturbance profiles closest to the optimal ones. 
Another option is to solve the receptivity problem for distributed generators upstream 
from the starting point, x in , to find their distribution, leading to the optimal disturbances. 
The next option is to design a disturbance generator that provides impact on the flow 
inside the boundary layer instead of a perturbation of the near-wall region only. For 
example, it might be a focused laser beam proceeding from the wall, where it could be 
delivered by a fiber-optic system. These fundamental issues should be addressed in 
future research programs on the application of bypass transition mechanisms to 
separation flow control at low Reynolds numbers. 
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Appendix 


For purpose of consistency of the formulation, we repeat the main features of the 
numerical scheme described in Ref. 5. Matrices A, Bg, Bp and B 0 in Eq. (7) are as 
follows: 





O' 


( 0 

0 

~P 

0^ 

0 

0 





u 

0 

0 

0 


-p 2 

~ u y 

0 

0 

V 

u 

0 

0 

> B 0 - 

0 

-2 Vy-P 1 

-pv 

0 

v0 

0 

u 

0; 


l 0 

0 

-Vy-P 2 




( 0 

-1 

0 

0 > 


^0 

0 

0 

o 6 

Bi = 

-V 

0 

0 

0 

, Bo = 

1 

0 

0 

0 

0 

—2V 

0 

-1 

z, 

0 

1 

0 

0 


l 0 

0 

-V 

0; 


^0 

0 

1 

0, 


The scheme utilized for Eq. (7) in the streamwise direction is 


(Af)‘ -(Af)° = Ax[(B 0 /)‘ +(B 1 y,)‘ +(B 2 


hAf)" +1 -2(Af)"+i(Af)" 1 =Ax 


(B 0 /r i + (B 1 /,,r 1 +(B2/vd 


\/?+l 


, n > 1 


where n stands for the step number along the coordinate x. At each streamwise position, 
the one-dimensional boundary-value problem is solved using a spectral collocation 
method based on Chebyshev polynomials. Usually, we used 100 intervals along the 
coordinate x and 100 Chebyshev polynomials for the solution approximation. 

To find the optimal perturbations corresponding to the maximum of the energy ratio 
in (9c), the forward solution of Eq. (7) is accompanied by the backward solution of the 
adjoint problem. 5 ’ 6 In the present work, the adjoint system was discretized, i.e., the 
discretized adjoint equations were employed. Another approach is to utilize the adjoint 
form of the discretized forward equations (see discussion in Ref. 5). An arrangement of 
the iterations is described elsewhere. 5 ' 6 Usually, 2 to 3 iterations were enough to achieve 
convergence at the 0.1 percent level. 
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